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quark impact factor is included in an appendix. 
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1. Introduction 

In this paper we report on the calculation of the strong radiative corrections to 
the process 

p + p — >b + b + X. (1.1) 

This reconsideration is prompted in part by the fact that measurements at the 
Fermilab Tevatron indicate that the b cross section lies above the predic- 

tions P, H from O(af) perturbation theory 0, [|. Similar excesses have been re- 
ported recently in the photoproduction b cross-section ||, Q, and in b production in 
77 collisions |H| [TTJ. Two remarks need to be made about this situation. First, the 
transformation of the observed experimental cross section to a 6-quark cross section 
requires the inclusion of fragmentation and decay corrections which are subject to 
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theoretical errors. Although these corrections may be smaller when we consider the 



6-jet cross section JT2| rather than the 6-quark cross section, experimental results for 



the 6-jet cross section still indicate an excess over theoretical predictions Sec- 
ond, the fixed order perturbation theory description of 6-quark production receives 
large corrections in order which call into question the validity of perturbation 
theory. We may simply conclude that, in this energy and mass range, the theory 
is not very predictive and ignore the disagreement with experiment as an inadequacy 
of the theory. A more constructive approach, which we adopt in this paper, is to try 
and resum the large corrections. 

The heavy quark cross section is calculable because the heavy quark mass M is 
larger than the QCD scale A. As first observed in ref. the large corrections in 0(a^) 
come both from the threshold region s ~ M 2 and from the the high energy region 
s ^> M 2 , where s is the normal Mandelstam variable for the parton sub-process. 

The resummation of the threshold soft gluon corrections is best carried out in 
a;-space, where u is the moment variable after Mellin transform with respect to s. 
In tu-space the threshold region corresponds to the limit uj — > oo and the structure 
of the threshold corrections is as follows 



a 



oo 2j 



=1 k=l 



Detailed studies |IJ| indicate that such threshold resummation is of limited impor- 
tance at the Tevatron where 6-quarks are normally produced far from threshold. 
Resummation effects lead to minor changes in the predicted cross section and only 
a slight reduction in the scale dependence of the results. 

At high energy, heavy quark production becomes a two scale problem, since 
S ^> M 2 ^> A 2 and the short distance cross section contains logarithms of s/M 2 . 
As we proceed to higher energies such terms can only become more important, so 
an investment in understanding these terms at the Tevatron will certainly bear fruit 
at the LHC At small u the dominant terms in the heavy quark production short 
distance cross section are of the form 

2 oo oo , 

*~^5>&e). (1.8) 

For j = we obtain the leading logarithm series (LLx). Resummation of these terms 
has received considerable theoretical attention in the context of the fc^-factorization 
formalism [[Uj-|2T|. Although some work quantitative nature has been performed p2 



23| , so far the connection with the precise structure function data has been indirect 
since the unintegrated gluon distributions are not measured directly. In addition, 
both the difficulty in reconciling /^-factorization with structure function data from 



HERA pi[-||29||, and the theoretical uncertainty in the face of large [RDJ destabilizing 
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sub leading corrections [31]— [34] have hampered this program. However significant 
progress in understanding the origin and resolution of these instabilities has been 
made recently f3"5(|-[3"5[|, and resummed parton densities in the small x region can be 



determined from HERA data [39|. Consistent resummed calculations of heavy quark 
production cross-sections are thus finally a possibility. 

The aim of this paper is to collect together the relevant results for a quantitative 
calculation of the total heavy quark production cross section. Our formulae include 
the known fixed order results, as well as the resummation of the leading tower of 
high energy logarithms. The approach to this resummation is as follows. The 0{a 2 s ) 
and 0(a%) total cross sections are well known at the parton level. Using the results 
of the /cT-factorization, valid at high energy, we can resum the logarithms of s/M 2 
by determining the leading coefficients of eq. (|1.3|) for k = 2, ... ,00, resumming 



the singularities using the techniques of ref . |38| , and then combining with resummed 



parton distributions |39| to give a high-energy improved short distance cross section. 



This paper deals exclusively with results for the total cross section. This quan- 
tity is not experimentally accessible with current detectors which typically have a px 
threshold below which 6-quarks cannot be observed. The technical problems of ex- 
tending this analysis to less inclusive measurements may be formidable. The present 
analysis is intended to give us an understanding of the size of these terms, before 
undertaking the extension to experimentally more realistic quantities. 

The structure of this paper is as follows. In section || we consider the descrip- 
tion of the direct component of bottom quark photoproduction both in fixed order 
perturbation theory and in the framework of /c^-factorization. We use the latter 
formalism as a technical device to resum the logarithms of s/M 2 and cast the final 
result in the form of a high-energy improved collinear factorization formula. Sec- 
tion [3] repeats similar steps for the hadroproduction of heavy flavours. In practice, 
the treatment of the photoproduction of heavy flavours is more complex than de- 
scribed in section ^ because of the hadronic (resolved) component of the photon. 
The resolved component of the photoproduction cross section should be treated by 
the methods of section |3|. Similar remarks apply to heavy quark production in 77 
collisions. Section |] presents some analytic estimates. In section |5] we draw some pre- 
liminary conclusions. A complete description of the calculation of the heavy flavour 
impact factors is presented in appendix [A]. 

2. Photoproduction of heavy flavours 
2.1 Fixed order perturbation theory 

In this section we report on the results for the direct component of the photopro- 



duction of heavy flavours in fixed order perturbation theory f4"0|| . We first define 



reduced cross sections S with the mass dimensions removed, both at the hadronic 
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and partonic level, 

S(p) = M 2 a{M 2 ,S) 



J 7J I Pi 



EL 
M 2 



M' 2 a 7j (s,M 2 ,ii 2 ^ 



(2.1) 



Here S and s are the squares of the total hadronic and partonic centre of mass 
energies, respectively, M the mass of the heavy quark, p = 4M 2 / S,p — AM 2 / s = p/x 
and p is the renormalization and factorization scale. We can further express the 
partonic cross section as 



Pi J^2 ) =a em a s (p )/ 7j (p,— ) " 



M 2 



(2.2) 



The photon-hadron cross section is then given by the collinear factorization formula 



E(p) = a em a s (p" 



1 da; 



a? 



x' M 2 



F j {x,p 2 



(2.3) 



where F j (x,p 2 ) = xq(x, p 2 ), xq(x, p 2 ), xg(x, p 2 ), j = q,q,g are the parton momen- 
tum densities, renormalized in MS at scale p. The coefficient functions f 7 j have a 
perturbative expansion which, following ref. [40|, we write as, 



fyj P 



LL 

M 2 



f^(p)+^a s (p 2 ) 



4- } (p) + in 



(!)/ 



2 1 



M 2 



+ 0(a 2 ). 



(2.4) 



The functions / 7J - depend on the charge of the quark which interacts with the photon. 
To make these dependences explicit we further define the quantities, 



As Pi 



EL 

M 2 

P 2 



e Q c 19 Pi 



p_ 

M 2 



M 2 



A<7 ( Pi "77? ) = A? ( Pi 7/77 



M 2 



CnC. 



IS 



( p ' w) + e « d ™ ( p ' m*) ' (2 ' 5) 



where e q is the charge of the incoming light quark, while cq is the charge of the heavy 
quark. Results for the perturbative expansions of the functions, c ig , c iq and d ig are 



given in ref. EC 



We can also define the moments f w of all functions /, 



dpp- L f(p) 



F^p' 



dx x^-'Fix^p 2 ) 



[2.6) 



Taking moments of eq. ( |2.3|) we get 



dpp" l H(p) = a em a s {p 2 ) ^ fw"(j^ F " 



(P 2 



(2.7) 
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The moments of c 79 ,c 7g and d ig are easily obtained from the results of 
expressions for the partonic cross-sections £ are then 



. The 



M 2 



c (i) +] n fiL)cW 



a em a^4vr ( e 2 . 



7g u 



75 w 



+ 



+ e 



Q 



c (i) + i n f ii_ ] cW 



(2.J 



The form of these results in the high energy limit are then readily derived 
uj -> 

a 



15 



: as 



V79 



JJ_ 

M 2 
M 2 



7 / 41 _ 7 ^ 
9 V 27 9 n M 2 / 



+ 



41 7 ^ 



~7re Q a fl a em 9^ M 2 J uj 



a 

- + 



where as usual a = N c as/ir, and V c = N 2 — 1. The d 19 terms are non sin; 
thus contain no uj = singularities. 



(2.9) 
det, and 



2.2 High energy behaviour 

We now turn to the high energy behaviour of the photoproduction cross sections 
which contain logarithms of s/M 2 and need to be resummed. This resummation is 
performed by using the k-factorized expression for the cross section 



_ j tC IT V CC ' 



(2.10) 



=9,1,9 

The functions JF- 7 are the unintegrated parton distribution functions and S 7 j are the 
(lowest-order and gauge-invariant) off-shell continuations of the parton cross sections. 
The diagrams necessary to calculate £ are illustrated in figure [I] where k 2 = xp 2 + k. 
We can simplify eq. (|2.10|) by taking the moments with respect to the longitudinal 
variable, cf. eq. Q2.6|) , to undo the convolution: 



d 2 k 



71 



Jl 



M 2 ' M 2 



^(k 2 ;p) 



(2.11) 



3=1,1,9 

j£ is the impact factor which expresses the response of the quark-antiquark pair 
as a function of the transverse momentum of the incoming light parton. A further 
simplification can be achieved by defining Mellin transforms with respect to the 
transverse momentum squared variables: 



d 2 k /k 2 



-oo 



7T 



d 2 k 

ri? Vm 2 



7 T (V_ f_ 

w \ M 2 ' M 2 



(2.12) 



5 





Figure 1: Diagrams for the direct component of heavy quark photoproduction at high 
energy. 

where the dependence of the Mellin transforms on /i 2 has now been suppressed. 
The Mellin transform of T is defined on the fundamental strip a < Re(7) < f3 
corresponding to the assumed behaviours 

^(kV) ~ (k 2 r-\ k 2 ^oo, 

The existence of the integrated parton distributions 

F^ 2 )= dk 2 ^(k 2 ^ 2 ) (2.13) 
Jo 

implies that (3 > 0. For the impact factor J the Mellin transform is defined in the 
strip < Re(7) < 1. The inverse Mellin transforms are defined as, 



J "{w'w)= h^iy-w) J " H - < 2 - 14) 

where C is a contour lying in the fundamental strips defined above. Inserting these 
expressions in the representation, eq. ( j2.11 ), and performing the transverse momen- 



tum integration we then have 

k=q,q,g JC V P 7 

To understand how this result may be used to resum high energy logarithms, 
consider first the unintegrated gluon density .^(k 2 , solution to the BFKL 

equation: the Mellin transform ^(7) then has the form 

wl7j 1 - >(7) ' 1 J 

where Xo(t) = 2^(1) — ^(7) — ^(1 — 7) is the Mellin transform of the lowest order 
BFKL kernel. Assuming that the asymptotic behaviour (k 2 ^> /i 2 ) is determined 
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perturbatively (rather than by the properties of J 7 ^), the rightmost singularity of 
^(7) is then a simple pole in the fundamental strip at 7 = j 8 (a/u)), where j s (a/u) 
is defined implicitly through the relation 

and the condition 7 S (0) = 0. When expanded perturbatively we have 

It follows that in the approximation in which we keep only the residue of this simple 
pole, the inverse Mellin transform, eq. fl2.14|) is given by 

In obtaining eq. ( |2.19| ) we have taken the limit k 2 > /x 2 so that we can neglect the 
contributions from non-perturbative singularities in the input distribution !F^\^), 
as well as higher twist perturbative contributions coming from solutions to eq. (|2.17|) 
to the left of j s (a/uj)). Substituting eq. Q2.19D into the definition eq. ( |2.13|) of the 
integrated gluon distribution we get 

Ft(f) = - n /.ff^'U ' h. (-)) ■ (2-20) 

High energy factorization makes a definite prediction for the evolution of the inte- 
grated gluon distribution in terms of the function 7 S : 

^(M 2 )=^J i^V)- (2.21) 

In a similar way, if the integral eq. ( |2.15| ) may also be dominated by the simple 
pole at 7 = j s (a/uj) (because the impact factor is regular in the fundamental strip), 
we find that the gluonic contribution to S w is 

T»(«/ w ) 



aw fM^~'- { " M 



J W5J)lj>J ^V) + -- (2.22) 

where in the second line we substituted eq. (|2.20|) for the integrated gluon distribu- 
tion, and in the last used the fact that at high energy (i.e. uj — > 0) we may write 

7^(7)=j(7)+OM, (2.23) 
since the impact factor ^(7) is free of singularities and terms of 0{uS) are sub leading. 
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This simple argument is however not sufficient to define the high energy cross- 
section, since if we redefine the unintegrated gluon distribution through a scheme 

change ^(7) -> ^(7)^(7), the function x(l,a s ) = Xo(l) + a s Xi(a s ) H only 

changes at NLLx [gl], 0: 

xd) ■ 



X(i) - a s b x{j)-^fou{l) + 0(a 2 s ) . 



(2.24) 



where, as usual, b = IL-^JUJ. . Consequently unless the BFKL equation is treated 
at NLLx, the effect of the impact factor can always be removed by a scheme change. 
Fortunately the NLLx kernel Xi{l) is now known and techniques to stabilise 
the expansion have now been developed which make meaningful NLLx resummations 
possible |PS||. All that remains is to specify the particular factorization scheme in 
which the impact factor is to be calculated. Impact factor calculations such as those 
described in appendix [A], based on the evaluation of off-shell amplitudes implicitly 
employ a "fc^-factorization scheme" . For the photoproduction of heavy quarks, the 
result of such a calculation is jn| [T^] , cf. eq. ( |A.32j ) 



J (7) 



7T 7 — 57 



5(1-7,1-7)5(1 + 7,1-7). 



(2.25) 



33-27 

We note parenthetically that the lowest order expression for ^(7) including the full 
uj dependence has been given in ref. ||19|| . In our notation it is 



J%{l) = e 2 Q a em a s 4"tt B{1 - 7 + u, 1 - 7 + + 7, 1 

[14 + 20c; + 9lu 2 + cu 3 - 7 (10 + 7u + cu 2 )] 
X 7(3 + u){2 + a;) (3 - 2 7 + 2u) 



7 + x 



For small 7 



J (7) ~ e 2 Q a em a 



7n 417T 



9 + 27 



7 + 0( 7 2 ) 



(2.26) 



(2.27) 



This can be related to the result in MS factorization through a scheme change ^(7) = 
5(7), where R^) is a universal function which fixes the normalization of the MS 
gluon distribution f43|. Explicitly 



r(i-7)x(7) 

(r(l + 7)(-7) X '( 7 )) 



lV2 



x exp 



7 ^(1) + 



dx 



(^(l)-^(l-x)) 



so for small 7 



1 + ^Cs7 3 



^C47 4 + fC57 5 + 0(7 6 ) 



where ( n is the Riemann zeta function. In MS factorization we thus have 

jJ 9 M=j(l)R{j) + 0{u), 



(2.28) 



(2.29) 



(2.30) 
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Note that once we have fixed the factorization scheme, the leading order calculation 
of the impact factor is sufficient for a consistent NLLx calculation: if we knew the 
impact factor at NLO, we would also need the function \ a ^ NNLLx for a fully 
consistent calculation. 

It remains to add in the quark contribution. At LLx 7 OT ~ r dggi where r c is 
defined as the ratio of colour charges: 



V r 



2W C A 9 



(2.31] 



when N c = 3. It follows that a quark may turn into a gluon provided one includes an 
additional factor of r c . However ^ qg is NLLx, so if the gluon turns back to a quark it 
costs an additional power of a s . Consequently at LLx the (bare) unintegrated quark 
and gluon distributions are related by 



^(k 2 , / , 2 )=r c (^(k 2 , / , 2 )-^)(k)). 



(2.32) 



It is then straightforward to prove, along the lines of similar arguments presented 
in ref. ||43|| , that at LLx the impact factor for quarks is related to that for gluons: 
specifically in MS 



7 J2(7) = 7^(7) = r c (j(l)R(j) - j(0)) + 0(u) 



(2.33) 



Putting all these results together, we find that at high energies the partonic 
cross-sections in MS factorization are given by 



V7S 

E 7 , 9 



J 17s 

S 7 5 = 



R[j, 



( M 2X ls ^ u ^ 



J 7 



R[i. 



M 2 
"7 



2\ 7s(a/u) 



i(o) 



(2.34) 



2.3 Double leading resummation 

Using the expansions eqs. (p7T8|) , (|2~27|) and ( ^29|) 



J 7 



~ e a em a s 



7n Alir a 



-) 

us) 



(2.35) 



It follows that the high energy results (|2.34j ) are consistent with the high energy 
expansions (|2.9|) of the fixed order results ( |2.8| ) in the region in which they overlap. 
Consequently we may construct "double leading" expansions of the high energy cross- 
sections by combining the two, provided we take care to subtract the terms (|2~9l ) to 
avoid double counting. For the case of the gluonic contribution to photoproduction 
such a matched calculation has been already presented in ref. T^| . 
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So overall the result for the high energy cross-section in double leading expan- 
sion is 



e Q a em a s 



+ 



J 17. 



UJ 



.(1) 



7- - 



— -j(Q) l + -ln— --j'(O) 



^79 



^7<? 



,a:47r 



"79 u> 



+ In 



M 2 
p? 

M 2 J 



M 2 

7s(»/w) 



c (1) 

75 ^ 



+ 



M 2 



a; 



"79 u> 



+ 



+ e n a em a s 4:7ra s 



+ r c 



J 7 



c (1) + In 



a 

UJ 

M 2 



M 2 



c (1) 



#7 



M 2 

7^ 



+ 

7s KM 



i(0) 1 + -In 
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;/(0) 



(2.36) 



where j(0) = ega em a s 77r/9 and j'(0) = e 2 Qa em a s AlTi /27 '. The final step is then to 
implement the NLLx resummation by replacing ~f s (a/oj) with a subtracted anoma- 
lous dimension, as described in |39|. There is then an ambiguity in the treatment 
of the double counting terms, just as in the resummed double leading expansion of 
the anomalous dimension: we can choose whether to resum the full series of sin- 
gularities eqs. ( 2.34j) , or whether instead to omit the double counting terms from 
the resummation, so as to leave the 0(a^) terms unchanged. This cannot be re- 
solved at NLLx, and will thus be useful as an estimate of the overall uncertainty in 
the procedure. 



3. Hadronic production of heavy flavours 
3.1 Fixed order perturbation theory 

Before considering high energy resummation we first describe the results in collinear 
factorization for the total cross section for the hadropro duct ion of heavy quarks of 
mass M. The hadronic cross section is factorized as 

M 2 a(M 2 ,s) = = W ^(x 2 , , (3.1) 

z — ' J xi x 2 \{xix 2 ) M Z J 

where the sum runs over all i, j = q,q,g, = £ji and p = AM 2 /S. Following 
ref. || the partonic cross section is given by 

M 2 a tJ { S) M 2 , /i 2 ) = % (p, jj^j = a 2 s (p 2 ) f i3 (p, , (3.2) 
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where p = 4M 2 /s = p/ixiX^), s being the partonic centre-of-mass energy. The 
functions = fji have a perturbative expansion 



fij P 



pS_ 

M 2 



flf{p)+^a s {p 2 



fS\p) + fS\p)^ ' 



M 2 



+ 0{a\ 



(3.3) 



The functions / have been calculated in perturbation theory to NLO in ref. ||, f(°\ 
p 1 ' exactly and as a numerical fit. Note that while T, gg and T, g g = begin 
at 0(a 2 ), f^gg = 1b gq - only start at 0{aP s ) and T, qq = S 55 (currently unknown) start 
at 0(af). Moments f u of all functions / may be defined as in eq. ( |2.6| ). Taking 
moments of eq. ( pM| ) we get 



= f dp p^Zip) = a 2 (p 2 ) £ f l3 u Fl(p 2 )F^) . 

% 3 



(3.4) 



With this normalization, we find the following results [ 175 1 in the high energy limit 
as uo — > 



(o) ~ 



(!) ~ 



'99 u 



V r V 15 c 18 N, 



O(uo) 
41 1 M 



f, 



99 ^ 



:N r - 



7 1 



2nV c \15 18 NJ uo 
and hence that (in the same notation as eq. ( |2.9|) ) 



+ 0(1) 



(3.5) 



^l 9 



v c 

2N 2 



_ T_l_ 

15 c 18iV r 



1-2- In At +2 , 

u M 2 uo \ 225 



a /154 41 1 



+ 



4 Ar 7 1 \ a , u 

— N c In h 

15 c 



18 N r J uo M 2 uo \ 225 



a A54 41 1 

c ~ 54 iv; 



+ 



(3.6) 



3.2 High energy behaviour 

At high energies the total cross section for the hadroproduction of a (heavy) quark- 
antiquark pair is dominated by the (Regge) gluon fusion process of figure |[ In the 
large centre of mass energy limit the perturbative expansion of the short-distance 



ki=ZiPi+k p 3 
-> 




ka=2" a p a -k' 

Figure 2: Diagrams for heavy quark hadroproduction 



11 



cross section contains large terms of order o? s [a s (M 2 ) log(s/M 2 )] n which need to be 
resummed to all orders. This resummation is performed by using the k-factorized 
expression for the cross section, which here takes the form 

S (P) = E / — ———^(k,k>,M 2 ,z 1 z 2 s)F(z 1 ,k 2 ^ 2 )x 

x P(z 2 ,k' 2 ,» 2 ), (3.7) 

in which Ti are the unintegrated parton distribution functions and S is the (lowest- 
order and gauge-invariant) off-shell continuation of the parton cross section. Taking 
Mellin moments with respect to the longitudinal momenta (see eqs. (|2.6| ) and 
again undoes the convolution: 



f (3.8) 



A further simplification is obtained by Mellin transform with respect to the trans- 
verse momenta: 

^ = Y.j c ^0 (jf) 3*M ^'(7i>7 2 ) , (3.9) 

where the Mellin transforms of the unintegrated parton distribution functions are 
defined as in eq. (|2.12 ) and the impact factor 

, r d 2 k d 2 k' ( k 2 \ 71 / k' 2 \ 72 TTlj ( k k' \ , , 



vrk 2 vrk' 2 \M 2 J \M 2 J u \M' M 

When we pick up the residues of the two perturbative poles eq. (|2.16|) in the 71, 72 
integrations, we find in a sequence of steps analogous to those leading to eq. ( |2.22| ) 

- E (% Bf^ (* © . * O) (f ) ^ ■ w 

ij 

where F^(fj, 2 ) are the integrated parton distribution functions at scale /z 2 . 

It thus remains to determine the impact factors 1/^(71,72). This requires the 
calculation of the moments of the heavy quark production cross section with two 
incoming off shell gluons. Details of this calculation may be found in the appendix. 
The answer in our normalization, and for u = 0, is 
2 

#0(71, 72) = -j-##(7i, 1 - 7i)5(72, 1 - 72) x 



x <UiV c 



B(3 - 71 - 72, 3 - 71 - 72) 
(1 - 71 - 72) 
B(3 - 71 - 72, 3 - 71 - 72) 



+ 



(1 - 71 - 7 2 ) 3 B 2 (l - 71, 1 - 72) 

2 r(2 - 7i)r(2 - 72 )r(2 - 7l - 72 ) (7 - 5(71 + 72) + 37172) 

N c r(4-2 7l )r(4- 2 72 ) (1-71 -72) 



(3.12) 
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This calculation has previously been reported in ^TJ. Note that eq. ( [3.12j ) is in 
disagreement with ref. |2(| eq. (2.10)] and ref. |21], eqs. (3.14)]. 

When uj = H^ji, -y 2 ) has a triple pole at 71 = 1 — 72, which splits into a simple 
and a double pole when u is non zero. In the vicinity of 71 = 72 = 1/2 one finds the 



singularity structure [IS 



7172^(71,72) = -T7- -r 
V c 6 



(1 - 71 - 72) (1 - 7i - 72 + 



(3.13) 



This singularity comes from the propagator for the s-channel gluon in the third 
diagram in figure ||: it is thus intrinsically non abelian, in the sense that it is not 
present in the high energy direct photoproduction cross-section given by the diagrams 
of figure |]. The singularity can drive a substantial rise in the hadroproduction cross- 
section at high energies, as we will show in detail in the next section. 

Since away from these singularities -£^(71,72) is regular in u, we may write 



7172^(71, 72) = h(ji, 72) + 0{u) . 
In the perturbative limit for small 71 and 72 we then have 



^(71,72) 
so that 



V 



4 

15' 



-Nr. 



7_l_ 
18 ~N f 



* f 154 
+ (7i+72) l^N c 



41 1 

54 ~N F 



+ 



h(Q,Q) 



k<\: 1 4 ^ 7 1 

15 c ~~ 18iVc 



7=0 



V e 

= h'(0,0) 



irat 



154 



V r V 225 



*r 41 1 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



As in the photoproduction case we must take care to identify the factoriza- 
tion scheme correctly. In particular to change eq. (|3.12j ) from the "/^-factorization 
scheme" to MS factorization now introduces a factor of R^i) R^) , with R the 
MS-gluon normalization factor eq. (|2.28|) , so that in MS 

7x72^(71, 72) = M7i, 72)^(71)^(72) + 0(u) . (3.18) 

Furthermore arguments similar to those in ref. [^3j may again be used to determine 
LLx high energy heavy quark production cross sections with incoming quarks in 
terms of those with incoming gluons: again in MS (compare eq. ( |2.33| )) 

7172^ (7i, 72) = r c (%i, l2 )R{ lx )R{ l2 ) - h( 7l , 0)R(n)) + 0(u) , 
7i72^ 9 (7i,72) = r c (/i(7i,72)i?(7i)^(72) - h( ll ,0)R( ll ) - 

- h(0, 72)^(72) + h(0, 0)) + 0{u) , (3.19) 



and similarly for the various other combinations. 
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Putting all this together, in the high energy limit the partonic cross sections for 
hadropro duct ion of heavy quarks are 



OJ 



10 



h is - ,7- - R z hs - K 



2\ 27 s (a/aj) 



^gq — ipgq 



<i \ / (/ \ \ „ 2 / Saw/ 



^ 7 



2\ 27 S (a/a;) 



-M,.e).o)«(-e(f 



a; / / \ fi- 

2~ 



_ yao _ y<?<? _ y<?<? 

J a; — — ^a; — ^a; 



2\ 27„ (a/a;) 



r/ \ / r/ \ \ ^ 2 / / (I \ \ ( M 2 



2/i 7 . 



a; 



,o) j r(t»( 



a\\ (M 



2\ 7s{a/w) 



+ 



+ h(0,0) 



(3.20) 



3.3 Double leading resummation 

The perturbative expansion of the MS partonic cross-sections eqs. ( |3.20| ) may be 
derived using eqs. flHSD , ( ggg ) and (|3~T5|) : 



h Is - ,7. F ))& (is - 



- /i(0, 0) + 2^/i'(0, 0) + O fa 2 , 

» (* O - °) fl 8) - "< -°> + >>°> + («• ■ < 3 - 21 > 

with h(0, 0) and h'(0, 0) given by eqs. ( |3.16|) and (|3.17|) , consistent with the fixed order 
results eqs. ( |3.6|) . Combining the fixed order cross-sections with the high energy cross- 
sections, taking care to subtract the terms eqs. (|3.6|) to avoid double counting, we get 



a 



+ 



s ] Jggu 



Is 



f (i) + i n ( ifl^ f(D 

# ( 7 . ( 



a\\ / M'' 



2 \ 27s (a/a;) 



+ r c 



fr(0, 0) f 1 + 2- In ^ - 2-/i'(0, 0) 
{tea, [/<JU In (£) } 



M-.l-).-.(-))'v J (-..(-l)l^ 



a; 



2\ 27 S (a/a;) 
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OJ 



10 



2\ 7a(o/w) 



. a , M 2 a , , . 
/i 0,0 - In — - -/i'(0,0) 



E« = E* = a 2 /<& + 4vra s 



f ( y +ln 

J qq u! * 111 



+ 



+ ri 



a; 



2\ 2-y s (a/uj) 



2\ 7s(a/a>) 



M%(-)-%(-))^(%(-)) 



a-' 



a; 



2 \ 27 3 (a/w) 



2\ 7»(<»/w) 



+ Mo,o, 



. (3.22) 



The NLLx resummation [^9] may be implemented just as in the photoproduction 
case, with similar ambiguities in the treatment of the double counting terms: here 
this ambiguity should be relatively unimportant at high energies since it has no effect 
on the singularities of h( i y s (a/u}), , y s (a/u)). 



4. Asymptotic behaviour at high energy 

The aim of this paper has been to calculate and collect results so that a matched 
numerical calculation of heavy quark cross sections is possible. The matching is 
necessary so that the calculation correctly includes the known fixed order results as 
well as the leading tower of logarithms present at high energy. 

In the following we give some analytic results for the high energy resummation 
alone; these are in no way meant to supplant the complete analysis alluded to above, 
but instead serve to indicate some features of the high energy resummation. First, we 
will obtain results for the asymptotic behaviour of the gluon distribution function at 
high energy, keeping the coupling fixed. We consider two alternative choices for the 
anomalous dimension: an anomalous dimension given simply by the singular term 
in lowest order perturbation theory (so that x ~ 1/7 an d has no minimum) and, for 
contrast, the lowest order BFKL anomalous dimension (so that the corresponding 
X-function has a minimum at 7 = 1/2). Second, we will explore the consequences 
of these two forms of the anomalous dimension for the high energy behaviour of the 
heavy quark photoproduction and hadroproduction cross sections, expressing the 
results in terms of the integrated gluon distribution function. For simplicity, only 
the gluon contributions will be considered, since the quark contributions are smaller. 
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Last, we will attempt a more quantitative estimate of the expected enhancement by 
letting the coupling run and choosing 7 = a(- — 1), so that x{l) = 1/(7 + a ) an d 
the momentum conservation constraint at uj — 1 is satisfied explicitly. This form is 



a good approximation to the fully resummed anomalous dimension [By], 38. 39 1 and 
furthermore gives a good account of scaling violations in HERA structure function 
data . 



4.1 Gluon distribution 

In x space, the integrated gluon distribution at the heavy quark mass scale, M, can 
be obtained by performing the inverse Mellin transform of eq. ( p. 21 ) with respect 
to u 



duj ,, ( M 



-X(IW^+^G)') GV >< (41) 

where L x = ln(l/x), I = lia^M 2 /^). Here G°(u) = F^{^) is a starting function at 
scale fio which we will assume is benign in terms of poles. For simplicity, we will work 
entirely in the "/c^-factorization scheme", so that the factors of R (see eq. (|2.28|) ) are 
absorbed into the gluon distribution. 

The asymptotic form of G(x, M 2 ) depends on the form of the singular anomalous 
dimension j s (a/uj). If we take the simple form j s (a/uj) = a/ui, the uj integration may 
be done by a saddle point approximation pi , with the result 



(4.2) 



where the argument of G° is the position of the saddle point. 

If, on the other hand, the anomalous dimension r y s {a/uj) is more like the BFKL 
form, it is convenient to make the usual change of variables uj = ax(j), (see eq. ( |2.17|) ), 



G(x, M 2 ) = / ^L(-ax'd)) exp Mt)^ + 70 G (ax( 7 )) • (4.3) 
Jc 2m 

If we now specialize to a x-function which has a minimum at one half 

1 „ n\ ( i x 2 



xh)^x[- 2 ) + - 2 x"[- 2 )[l-- 2 ) +■■•, (4.4) 

so that, as L x increases, the saddle point in the 7 integration moves to 1/2 from 
below. Asymptotically we then find when L x 3> / 

IL~ 3/2 ( . (\\ 1\ „„ / fl 
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Now consider the direct component of the photoproduction cross-section in the 
high energy limit. From eq. (|2.22|) we have that 

E <" = 1 (isj exp ( aL + -» O (->• O) G ° (w) • (46) 

where L = lnl/p. The function ,7(7), defined in eq. Q2.25Q is free from poles in the 
interval < 7 < 1 and is therefore a smooth function in the region of the saddle 
points of the integrand. It follows that at high energy for the "perturbative" gluon, 
eq. §3), 

~ G(p,M 2 ), (4.7) 

where uJ = \JaljL is the position of the saddle point. This expression will break 
down when the argument of j approaches one, since j has a simple pole there, but 
should still be good around a/uJ ~ 1/2, where it leads to a small enhancement [|i~8[| . 
For the BFKL gluon, eq. ( |4.5| ), we find at high energy 

m^j(l)G(p,M 2 ). (4.8) 



2 t 

In both cases the rise in the gluon drives the rise in the cross-section, with only a 
modest enhancement over the leading order result, j(0)G(p, M 2 ). 

4.2 Hadroproduction 

We now consider the high energy behaviour of the cross section for the hadroproduc- 
tion of a pair of quarks of mass M. From eq. ( [3.11| ) we have 

^(-^(7)) exp(ax(7)L + 2 7 0M7,7) [G°(a X ( 7 ))] 2 , (4.9) 

where ^(7,7) = 7 2 if w (7,7), L = lnl/p and in the last line we again made the 
change of variables uj = ax(j) ■ In contrast to the direct photoproduction cross- 
section considered above, the asymptotic behaviour at large L will now be dominated 
by the singularities of ^(71,72) (see eq. ( |3.13| )), 



(1 -2 7 )(1 -2 7 + w) ; 
A 



IV 2 



.(1-2 T ) (l-27 + w) (l-2 7 + cu) 2 _ 



(4.10) 



with A = ira 2 -^-. 

s OVc 
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First consider the "perturbative" anomalous dimension 7 = a/ou. Because the be- 
haviour in uj is more intuitive we choose to work with the uj representation, eq. 
With this choice for 7, eq. ( }4.10| ) becomes 



, f a a 

K[-,~ 

\UJ UJ 



A 



UJ + Co>_ 
+ 



LU + LU-(UJ — UJ ) (UJ + — LU-) 2 { (uj — UJ + ) 2 (uj — CO-)' 
{ljj 2 _ — 3dJ_Co> + ) (3co>_CJ + — toV) 



(4.11) 



{uj + — Co>_) 3 [ lo+(uj — UJ + ) LU-(iU — UjJ) 

where 

[-1 ± VTT8d~] . . 

u = 2a, uj± = . (4.12) 

The dominant behaviour at large L and I depends on the position of the saddle point, 
ZJ = yj2al/ L relative to the rightmost pole, uj = oj . If the saddle point is to the 
right of the pole, uj > 2a, (which means 2aL < I), we integrate along the contour 
uj = ZJ + iz to obtain the result 

S(p) ~ — exp (ujL + 2=A rC°(cJ)l 2 / dz f(U + iz, uj , uj + , ujJ) x 
2tt V uj J 

x exp , (4.13) 



where / is given in eq. (|4.11|) and the quantity c, 

rr 3 \ 1/2 

k) (414) 

sets the scale of the gaussian fluctuations which determines whether the poles at 
ojq,u + and uj- should be considered close. 

Using the saddle point method in the vicinity of poles we need to consider inte- 
grals of the general form 

= <~£«fa^«p(-£), (4.15) 

with y = b/c and b, c 7^ 0: K n (—y) = (—) n K n (y). An evaluation of the first three of 
these integrals Ki(y), K 2 (y) and K 3 (y) may be found in the appendix ||. In terms 
of these functions the result is 

Ac 

S(p) = j^yjt X(uj,uj ,uj + ,uj_) exp(2v / 2^L) [G°(uj)] 2 , (4.16) 

where 

_ 1 f°° _ ( z 2 

X(uj,ujq,uj + ,ujJ) = — -= / dz f(uj + iz,uj ,uj + ,ujS) exp — - 

Cy/TT J_ 00 \ C 
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C\/TC 



u - u ) 



UJ+UJ- \ c 

(ZJ — u+) \ u)-{oj_ — 3toL 



c J lo + (lu+ — C<J_) 3 



+ K x 
-K 2 



C J c(u + — CO-)' 

(u; — U-) \ l) + (oj + — ?>ojJ) 



C I U)-(u) + — CO-) 3 



,uj — uj- \ uj- 



C J c(iU + — LU-) 2 

Expressed in terms of the gluon distribution eq. (|4.2f ), we find 



(4.17) 



- ^T^yi/T X(cJ, u , w+l w _) [G(VP, M 2 )] , (4.18) 



where we have used the fact that if x = y/p, L x = |L and thus 2\/2alL = A\/2alL x . 
As To — > co>o from above, the enhancement factor X(u7, c^o, u;_) grows until it 
reaches a maximum value at the pole. 

When uj passes the pole, we have to add the residue of the pole to the saddle 
contribution. The result for the pole contribution for 2aL > I is 

S(p) ~ ^exp(2aL + /) [G°(cu )] 2 

~ A^exp(2aL-/) [G(e^ /4a , M 2 )] 2 . (4.19) 

At high energies this contribution is dominant, since it grows as e 2aL = p~ 2a . Taken at 
face value, eq. (|4. 19| ) thus indicates a dramatic growth of the hadropro duct ion cross- 
section with increasing energy, driven not by the growth of the gluon distributions 
but rather by the perturbative singularity in the partonic cross-section. 

If instead we consider the BFKL anomalous dimension, the asymptotic behaviour 
of the hadroproduction cross-section is best determined using the 7-representation, 
eq. ( fO| ) and h w as given by eq. ( f4.1Up , The function ec L- ([1-40, has a minimum 
at 7 = 1/2, and S is thus approximately given by 

aA X »(l/2) f dl exp (a X ( 7 )L + 2 7 /) (4 2Q) 



o 2<Ki (7 - I - \a X (l/2)f 
The integrand now has a saddle point at 



and the fluctuations about the saddle are determined by c 7 = a/2/ (ax"L). The 
saddle point is always to the left of 1/2 which it approaches asymptotically as L 
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increases. The double pole is at a fixed distance to the right of 1/2. The behaviour 
at large L is thus given by 



S(P) 



ir ~~ ^ J eX P 

lb7rc 7 \ 2c 7 / 



: exp ax 



2a[x(l/2)]V27raL 
which when expressed in terms of the BFKL gluon eq. 



(f | ax 

~5[) becomes 



(4.22) 
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/ a X "(l/2) 




2 L 5/2 




K x(i/2) J 


Z 2 



(4.23) 



Clearly in the BFKL case, although there is some enhancement JT8[ , it is rather less 
dramatic than in the previous case. This is because when x(t) nas a minimum, the 
saddle point cannot travel past the pole at 7 = 1/2. 

In summary, because of the singularity structure of eq. ( |3.13| ) substantial en- 
hancements of the hadroproduction cross-section are possible at high energy, even 
when the gluon distribution function does not have the singular behaviour predicted 
by lowest order BFKL. 

4.3 A more realistic model 

In order to attempt a quantitative estimate of the enhancement in the hadropro- 
duction cross-section at the Tevatron and LHC, we take a simple model for the 
anomalous dimension function which is close to the one-loop value and satisfies mo- 
mentum conservation, 

1 



7(w) 



1 



(4.24) 



We also include running coupling effects. In a theory with a running coupling the 
product la is replaced by 



CM 



(4.25) 



If we take the leading order parameter A = 200 MeV (corresponding to the leading 
order coupling at the Z-mass a s (M z ) = 0.123), we find that a = 0.224. Further, 
taking // = 1 GeV from eq. (|4.25|) we find ((M) ~ 1. We shall adopt this value in 
our numerical estimates below, although the details will change if we use different 
values for /x and A. 

The modified anomalous dimension, eq. ( 4.24 ), only changes eq. ( |4.2j ) in a mi- 
nor way 



G{x^ 2 ) 



1/4 



(CM) 

(4tt)V2lJ 



3/4 



cxp 



G° 



CM 



(4.26) 
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Tevatron 


LHC 


L = ln(,S/4M 2 ) 


10.6 


14.5 


C(M) 


1.0 


1.0 


cJ=(2C(M)/L) 1/2 


0.43 


0.37 


uj = 2a/ (1 + 2a) 


0.31 


0.31 


uo + = -§- a + ((§ + a) 2 + 2a) 1/2 


0.26 


0.26 


uo_ = -\-a-{{\ + a) 2 + 2af 2 


-1.71 


-1.71 


c= [oJ 3 /2C(M)] 1 / 2 


0.20 


0.16 


K^o -5?)/c| 


0.61 


0.39 


|(o; + -aJ)/c| 


0.85 


0.68 


l(w--w)/c| 


10.6 


13.0 


X(u, u ,u+,u-) 


1.26 


2.42 



Table 1: Values of parameters at Tevatron and LHC, taking M = 5GeV, a(M 2 ) = 0.224 



More interesting is the effect on the hadroproduction cross-section. h u still has the 
form given by eq. ( |4.11|) but the positions of the single and double poles are now 
given by 

2a 



(4.27) 




*=(^Y /2 . (4.28) 



The saddle point is now at 

- f 2 < 
iii = I — 

L 

Some plausible numerical values of the parameters are given in table ffl. The leading 
pole is at u = u Q , closely followed by the double pole u + , while u_ is in practice of 
little consequence. 

We see from the table that To > ujq at both the Tevatron and LHC, so we will 
concentrate on this case. The more dramatic behaviour, eq. ( }4.19|) only seems to 
set in at yet higher energies, of the order of a few hundred TeV. The scale of the 
gaussian fluctuations which determines the closeness of the poles is now set by 

/ _ 3 \ 1/2 

c= (Mm)) 

The saddle point result is (using eqs. fl4.11|) and ( [4.17|) ) 

/ _3 \ 1/2 

S(p) ^ A \^^) exp (2V2C(M)L - 2C(M)) X(ZJ, u , u+, u;_) x 

x [G°(uJ)} 2 . (4.30) 



(4.29) 
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We can express eq. (|4.30|) in terms of the gluon distribution function, eq. Q4.26Q 
and obtain 

s(p) - a WmW x{p ' Wb ' ^ w-) [G(v ^' M2)] ' (4 ' 31) 

In order to estimate a X-factor we calculate the LO cross section with 

7172^(71,72) — ► h(0,Q) (4.32) 



where (cf. eq. (|3.6|)) /i(0,0) is given by eq. (|3.16|) . The result is 

-jt-1/2 r 3/4 

= M0,0) ^ [GtvfoM 2 )] 2 , (4.33) 

and hence the iT-factor is 

where A/h(0, 0) = 135/181 ~ 0.75. With the parameters in table [l], the numerical 
value of K is thus around 1.4 (2) at the Tevatron (LHC). These numbers indicate 
that resummation can give rise to sizeable enhancements. A full numerical analysis 
is required to show whether they actually do. 



5. Conclusion 

We have presented a matched formalism which may be used to provide accurate 
estimates of the heavy quark cross section at high energy. Simple analytic calcula- 
tions suggest that the high energy resummation can give large effects, even with a 
gluon distribution which has the form dictated by lowest order perturbation theory. 
The cause of this enhancement is an s-channel gluon which leads to the singular 
behaviour of the non-abelian contribution to the impact factor in eq. ( |3.13| ). The 
effect is therefore universal, in the sense that at high energies it will provide the 
leading contribution to hadroproduction, resolved photoproduction and even double 
resolved 77 heavy quark production. The persistence of the singular behaviour in yet 
higher orders remains an open question. Furthermore, detailed numerical work will 
be required to establish the quantitative significance of the enhancement at existing 
and future colliders. 
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Note added in proof. After the appearance of this paper, we were informed 



by Marcello Ciafaloni that he is aware that some of the equations in refs. pO , [21 



are indeed incorrect. An erratum correcting some of the equations in ref. |2lJ has 
now been submitted to Nuclear Physics B in May 2001. The result of Camici and 
Ciafaloni for the hadronic impact factor as corrected by their erratum to ref. [2T| is 
now in agreement with our eq. (|3.12|). 



We would like to thank Marcello Ciafaloni for correspondence on this matter. 
A. The hadronic impact factor 

We wish to calculate the off-shell cross section (impact factor) for the production of 
a heavy quark pair by a gluon fusion process, 

g(k 1 )+g(k 2 )^Q( Pi )+Q{p 3 ). (A.l) 

The incoming and outgoing momenta are specified with respect to lightlike vectors 
Pi and p2- We perform the Sudakov expansion of all four momenta, where k, k' and 
A are euclidean vectors in the plane transverse to pi and p 2 , 

h = Z\V\ + k , 
h = Z2P2 - k' , 

p 3 = (1 - x 1 )z 1 p 1 + X2Z2P2 + k - A , 

Pi = xizipi + (1 - x 2 )z 2 p2 — k' + A . (A. 2) 

The specific form of the incoming momenta, ki and &2 in eq. ( A.2|) gives the dominant 



contribution at high energy and leads to a gauge invariant cross section [|17], |T? 
The matrix element squared, \A\ 2 is calculated taking the polarization sums of the 
incoming gluons to be, 

k^k" k'^k /v 
£ £ A(^r (h) = -^r , £ £ x(h)eT(k2) = -j^- . (A.3) 

A A 

The dimensionless partonic hard cross section S is obtained from the matrix 
element, A, by integrating over the phase space, including the flux factor, {2u = 
4ziz 2 pi ■ P2), and multiplying by the square of the heavy quark mass, M 2 , 

M 2 r 

£(k,k',MV) = — J d^\A(k 1 ,k 2 ,p 3 ,p i )\ 2 . (A.4) 

Our aim in this appendix is to calculate the arbitrary moments of this cross-section 
with respect to the transverse momenta k 2 and k' 2 and the zeroth moment with 
respect to the variable v. We define the moments of the cross section as 

r°° d 2 k f°° d 2 k' ( k 2 y 1 /k /2 \ 72 

#471,72)=/ -rp; \t^\ 77? x 



vrk 2 J nk' 2 \M 2 J \M 2 J 

x r±(^Y^ k , M ^ y (A . 5) 



v \ v 
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The results of this calculation for uj = have been been reported in refs. [20, 21]. We 
are in disagreement with the final results reported in these references. 1 We follow 
the general calculational procedure of ref. |21[] and present this appendix because we 
believe that a complete description of the correct calculation is a useful addition to 
the literature. 

We start with \A\ 2 , the (off-shell) squared matrix element for the gluon fusion QQ 
production process given in eq. ( |A.1| ) calculated with polarization vectors, eq. ( |A.3| ). 
In an explicitly Lorentz invariant notation we have, 



\A(ki,k 2 ,p 3 ,p 4 )\ 2 

j^{ab) 



9 



2N, 



2V r 



(A.6) 



(pm) 2 



-i 2 



(pm)' 



(t — M 2 )(u — M 2 ) 



h 2 h 2 

(pm) 2 



J^nab) 



(pm) + 2 

2 



(pm)(p2Pi) , n (pm)(P2P3) 



+ 



(t - M 2 ) 

4(pm) 2 
l 



1 



M 2 ) 



(t- M 2 )(u- M 2 ) 

1 N ( (pm) _ 



{t-M 2 ) (u — M 2 ) J \ (pik 2 ) 

(pm) 



(A.7) 



(P2P4) \ 

mki)) 



+ 



21,2 
2 



s (pik 2 )(p 2 k 1 

(pm) + 4 



(pm)(p2Pi) _ A 

(t - M 2 ) s 



x 



/ x . A (vm)(v2Pz) , a 
(pm) + + — 

\U — M Z ) S 



where 

s = (ki 
and we have defined 



k 2 ? 



t=(ki~Pzf 



u 



A = 2{p lP2 ) 



(piPi)(p 2 pz) _ 2 (pm)(p2P4) 
(pm) (pm) 

+ k 2 + (p-iki) - (p 3 k 2 ) 



(h-p 4 f 



k 2(pm) 



(A. 



(A.9) 



(Pih 



+ 



2 (P2P3) 

** (P2h) 



(A.10) 



eqs. ( |A.6| )-( [A.1(J ) are identical to the expression of Catani, Ciafaloni and Haut- 
mann |18| except that we have performed the interchange p^ *-> p^, necessary in our 
notation, eq. (|A.1| ). It is also equivalent to the expression of Collins and Ellis fl~7|j . 



Explicitly, we disagree with eq. (2.10) of ref. §(| and eqs. (3.14a), (3.14b) of ref. Q. In addition, 
we believe that the description of the calculation in ref. |2lJ] contains errors or typographical mistakes 
in the final results cited above, in eqs. (3.3), (3.7), (3.10a), (3.10b) and in the appendix eqs. (A. lb), 
(A.lc), (A.5), (A.6), (A.9), (A.10) and (A.16). 
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Using the Sudakov parametrization of the final state momenta, eq. ( |A.2|) , we can 
write the two-body phase space of the produced QQ pair as follows, 

d$( 2 ) = -H- dxi dx 2 d 2 A 5(x 2 (l - xx)v - (k - A) 2 - M 2 ) x 

87T 2 

x5(xi(l -x 2 )v- (-k' + A) 2 -M 2 ) (A.ll) 

1 ,2X A 2 + M 2 2 \ 

- cTA5 z/ 7-Q ) ( A - 12 ) 



1 dx 2 l2Z ./ A 2 + M 2 



where 



-^^-^^'("-^j-^)- (A ' 13) 

q = k - k' , v = z 1 z 2 (2p 1 ■ p 2 ) = s + q 2 , 

A = A - kxi - k'(l -xi) , A = A - k'x 2 - k(l - x 2 ) . (A.14) 

It proves convenient to use both parametrizations of the phase-space as given in 
eqs. ( A.12Q and ( A. 131) . Different terms in the integrand are more easily integrated 
using one or the other forms of the phase space. For the phase space parameterized 
as in eq. ( |A.12j ) we find that 

A 2 + M 2 



X\(l — Xi) ' 

M 2 —t — xi(s + k 2 ) + (1 - X!)k' 2 + 2A ■ k', 
M 2 — u = (1 — xi){s + k 2 ) + xik' 2 - 2A ■ k', 
[(1 -xi)q- A] 2 + M 2 



x 2 



1 -x x )v 



x _ [x iq+ A]' + ^ 

X\V 

whereas for the phase space as described by eq. ( JA.13Q we have 

A 2 + M 2 

x 2 (\ - x 2 y 

M 2 — t = x 2 (s + k' 2 ) + (1 - x 2 )k 2 + 2A ■ k, 
M 2 — u = (1 — x 2 )(s + k' 2 ) + x 2 k 2 - 2A ■ k, 
[(1 -x 2 )q + A] 2 + M 2 



(1 - x 2 )v 



1 



x 2 q — A + M 2 /a \ 

1- Xl = [ — — . A.16 

x 2 v 

To simplify the formulae we define a reduced matrix element squared, D 

4 2 

\A(h,k 2 , P3 ,p,)\ 2 = ^ r D(k,k',A,x 1 ,M 2 ). (A.17) 
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Expressed in terms of the reduced matrix element, D, and after the inclusion of the 
other factors from eq. (|A.4j) , the Mellin transform of impact factor, eq. (|A.5|) , takes 
the form 

d 2 k d 2 k! 



-ZL (71,72) 



a 2 s M 2 
2K 



, MM 2 

du\ 

V v 



vrk 2 vrk' 2 



dx\ dx 2 d A x 



xuS(x 2 (l — X\)v — (k 
/ k 2 \ 71 / k' 2 x 72 



x 



- 



\M 2 J 



M 2 



-A) 2 - M 2 )5(x 1 (l-x 2 )u - (-k'+ A) : 
£>(k,k', A,x l7 M 2 ) . 



M 2 )x 



(A.18) 



We now perform a further separation of the reduced matrix element, eq. ( A.17|) , 
to isolate the sub-leading and leading colour pieces, corresponding to the abelian and 
non-abelian contributions respectively, 



D 



pa) + D m + D 



(3)^ 



(A.19) 



1 



(M 2 — u)(M 2 



D {2) = N r 



t) 

1 



M 2 -t 
2{B-C) 



(B + Cf 

k 2 k /2 

(1 - Xt - X 2 ) 



k 2 k' 2 S 



D {3) = N r 



2 

us 



M 2 - u 
((l-x 2 )k 2 + (l-x 1 )k' 2 

2 ((l-x 2 )k 2 + (l-x 1 )k' 2 - 



B 2 + C 2 
~ k 2 k' 2 

- k • k') 
k • k'V 



k 2 k /2 



The quantities B and C are given by 
1 (1 - Xi){l — x 2 )v 



B 

C 



2 
1 
2 



M 2 — u 
X\X 2 U u(\ 



+ 



v(\ — X\ 



x 2 



2s 



+ 



A • (k + k') 



x x -x 2 ) A • (k + k') 



(A.20) 
(A.21) 

M 2 -t 2, s ' (A ' 22) 

and the Mandelstam variables of the hard sub-process are given by eq. ( A.9|) . The 
non-abelian part has been further split into two terms, and D^ 3 \ each of which 
leads to a finite A integration. The integration over A is finite because of the gauge 
cancellations between the different terms in B and C which ensures that B,C — > A~ 2 
at large A. Note that eqs. ( ]A.20p are finite as k 2 and k' 2 tend to zero and that the 
singularities at small k 2 and k' 2 are only apparent. This can be established by putting 
the expression over a common denominator and using the identity (cf. eq. (TATTTD ) 



x 2 v 



x\v + k - k - 2 k • A + 2 k' ■ A . 



(A.23) 



The first step in the evaluation of the u = limit of eq. ( A.18 ) is to use the 
5-function in eq. (|A.12| ) or ( |A.13| ) to eliminate the variable v. The choice of which 
parametrization of the phase space to use is made so that the subsequent integration 
over A contains at most one non-trivial azimuthal integration ETI . 
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After combining the denominators using Feynman parameters and shifting the 
variable of integration in the normal way, the integration over A can be performed. 
We next perform the k 2 and k| integrations, which can all be evaluated using the 
following identities, 



00 d 2 k fk 2 \ 



- 



- 



yrk 2 \M 2 J (M 2 + £k 2 )' 
d 2 k x d 2 k 2 ( k 2 



(M 2 )-T 7 %«-7), 



(A.24) 



00 t2 



vrk 2 vrk 2 \M 2 



71 



k 2 
K 2 



12 



1 



00 t2 



d 2 k ± d 2 k 2 

7rk 2 ixk\ 
= (M 



M 2 



M 2 

(M 2 )-"^r 71 G 

71 / \ 72 

M 2 



(M 2 + eik 2 + e 2 k 2 



72 



r(7i)r(72)r(w - 71 - 72) 
r(n) 



(A.25) 



"^-71-72 



[M 2 + e(k!^^ 2 j 

r(n - 71 - 72)r(i - 71 - 72)r(7i)r(7 2 ) 
r(n)r(i- 7 i)r(i-72) 



(A.26) 



Finally we perform the integration over the Feynman parameter and over the variable 
X\ or x 2 - 

In the uj = limit, the result for the abelian piece is 



# {1) (7i,72) 



2 2 



x 



r v r(i- 7 i-72)r(7i)r(7 2 )x 

-S(l- 7l ,l- 7l )S(l- 72 ,l- 72 ) 



5(1-7!, 2 -71)5(1-72,2-72) 



na 



2 2 



(3-2 7l )(3-2 72 ) 

r(i - 71 - 72)r(7i)r( 72 ) 



(l + (l-7i)(l-7 2 )) 



71 r(l - 7 i)r(l - 72 ) 



V c N c 

x 471+72 ( 7 _ 5 ( 7l + 72 ) + 3 7l72 ) 
2 



64r(|- 7l )r(|- 72 ) 



- 2 2 



5( 7i , 1 - 7i)5( 7 2, 1 - 72 )r(2 - 7l - 72 ) x 



x 



T(2 - 7l )r(2 - 72) (7 - 5( 7l + 72) + 3 7l72 ) 



(A.27) 



r(4-2 7l )r(4- 2 72 ) (i- 7l _ 72 ) 

where we have written the same expression in three ways in order to facilitate com- 



parison with the wrong expression of ref. |HJ] and the correct expression of ref. QlS 



The last form is the form quoted in this paper. In the uo = limit, the result for the 
non-abelian pieces is, 



# (2) (7i, 



72 J 



-'<: . ~ 5(3-71-72,3-71-72) 



5 (3) (7i,7 2 ) 



not] 



'-AN, 



■4iV c 



(1 - 71 - 72) 

5(3 - 71- 72, 3 - 71- 72) r( 7 i)r(7 2 )r 2 (2 - 71-72) 



(1 - 71 - 7 2 ) ; 



5( 7 i,l-7i)5( 72 ,l-7 2 ), (A.28) 

(A.29) 



r(i- 7 i)r(i- 72 ) 
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The final result for H is given by the sum of all three contributions 

#o(7i,72) = # (1) (7i,7 2 ) + # (2) (7i,7 2 ) + # (3) (7i,7 2 ) 



(A.30) 



The gluon on-shell limit can be investigated by examining the residues of the poles 
when 7 — > 0. We define the function h 



h(ji, 72) = 7i72 #0(71, 72) ■ 
The result when we take one leg on shell is known, 72 — » (see ref. [15]) 



(A.31) 



not 

Vc. 



B(l + 71, 1 - 7i)^(l - 7i, 1 - 7i) x 



x 



Nr. 



2(2 - 7l ) 



1 (7-5 7 i 



: (3-2 7l )(5-2 7l ) 2iV c 3(3-2 7l )J " (A ' 32) 

Note that the colour suppressed piece of this expression gives the contribution in an 
abelian theory. It is thus related to the impact factor for the direct photoproduction 
of heavy quarks: in the notation of eq. (|2.23|) , j( 7 ) = -e 2 Q ^ 2V C N C /i (1) (7, 0). We 
can also check the zeroth moment of the normal cross section when both gluons are 
on shell, cf. eq. (|3.5|) 



h(Q,Q) 



■no, 

Vc. 



4 



1 7 



^15 iV r 18 



(A.33) 



As a final check we have evaluated eq. ( |A.18| ) for H numerically and find good 
agreement with eq. (|3.12| ) in the range < 71,72 < 1/2. 

The triple pole in eq. (|A.29[) may be of some phenomenological importance. It is 
easy to evaluate for arbitrary lo leading to a contribution, which for 71, 72 ~ 1/2 
may be written as, 

^(71,72) ~ ^ —r, \7~j (A.34) 

V c 6 7172(1 - 71 - 72j(l - 7i - 72 + u) 1 

B. The K functions 

Using the saddle point method in the vicinity of a simple pole we have to evaluate 
an integral of the form 



KM 



dz 



(b 2 + z 1 



cxp 



vrexp(?/ 2 ) [1 - Er%)] = vr exp(y 2 )Erfc(?/) 



(B.l) 



with y = b/c and b, c > 0. Similiarly, for integrals in the neighbourhood of a double 
pole it is useful to define 



K 2 (y) 



dz- 



b 2 -z 2 



(b 2 + z 2 ^ 2 
2 v ^-2yK 1 (y) 



ex P — 9 



(B-2) 
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12 3 4 

y 

Figure 3: Plot of the functions Ki(y), K 2 (y) an d K 3 (y). 
and near a triple pole 

= K 1 (y)-yK 2 (y). (B.3) 



The functions K 1: K 2 , K 3 have the following expansions: for small y 



while for large y 



l-^ + y 2 + 0(y 3 ) 



Ki{y) = 7T 



K 2 (y) = 2^[l-^y + 2y 2 + 0(y 3 )] 



K 3 (y) = 7T 



KM 

K 2 (y) 
KM 



1 - % + zy 2 + otf) 

7T 



y 



y 



I ! — + O ' 1 



2y 2 
3 



V 



(B.4) 



(B.5) 



A plot of the functions Ki, K 2 and K 3 is shown in figure 
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